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Abstract: Primordial perturbations with wavelengths greater than the observable uni- 
verse shift the effective background fìelds in our observable patch from their global averages 
over the inflating space. This leads to a landscape picture where the properties of our ob- 
servable patch depend on its location and may significantly differ from the expectation val- 
ues predicted by the underlying fundamental inflationary model. We show that if multiple 
fields are present during inflation, this may happen even if our horizon exit would be pre- 
ceded by only a few e-foldings of inflation. Non-Gaussian statistics are especially affected: 
for example models of locai non-Gaussianity predicting I/^lI ^ 10 over the entire inflating 
volume can have a probability up to a few tens of percent to generate a non-detectable 
bispectrum in our observable patch I/nl' I ^ 10. In this work we establish systematic 
connections between the observable locai properties of primordial perturbations and the 
global properties of the inflating space which reflect the underlying high energy physics. 
We study in detail the implications of both a detection and non-detection of primordial 
non-Gaussianity by Planck, and discover novel ways of characterising the naturalness of 
different observational configurations. 
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1. Introduction 

Inflation is the leading scenario explaining the generation of primordial perturbations. 
The imminent release of Planck data will provide cruciai new Information of the non- 
Gaussian statistics of primordial perturbations. This will help to efficiently discriminate 
between different realisations of inflation, potentially ruling out whole classes of models. 
Future Large Scale Structure surveys like Euclid will further improve our knowledge of 
the statistics of primordial perturbations opening striking new insights to the high energy 
physics behind inflation. 

The observable primordial perturbations probe physics during the last A'obs. ~ 60 
e-foldings of inflation after the horizon exit of the largest observable modes. Inflation 
may however bave lasted much longer, so that the observable part of the universe would 
constitute only a fraction of the entire inflating patch. Long-wavelength perturbations 
generated before the horizon exit of largest observable scales average to constants over our 
observable patch, effectively shifting the locai background field values. In adiabatic single- 
field models it is well known that the long-wavelength contributions amount to simply 
shifting the locai time coordinate |||]. The situation is different if more than one light 
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dynamical scalar is present during inflation, as in many models generating detectable non- 
Gaussianity. In this case the evolution is non-adiabatic and the long-wavelength modes in 
general shift both the time coordinate and the fields parameterizing isocurvature directions. 
Shifts in the isocurvature directions correspond to non-trivial changes of initial conditions 
which may significantly affect the observational signatures and especially the non-Gaussian 
statistics, irrespective of whether the isocurvature perturbations persist until today. 

The long-wavelength contributions are difFerent in different parts of the inflating region 
which generates variances in the physical properties of patches smaller than the entire 
inflating region. This leads to a landscape picture where the statistics of perturbations 
within a horizon patch depend on the location of the patch itself [||, as previously discussed 
in the curvaton scenario The observable perturbations can of course be described 

by an effective model covering only the last A'obs. e-foldings of inflation in our horizon 
patch. The background field values in the effective description are modulated by a random 
long-wavelength contribution specific to our location in the entire inflating patch. In the 
presence of isocurvature directions the effective model can be non-trivially related to the 
underlying fundamental inflationary model which covers the entire inflationary epoch A'tot. • 
Given a fundamental model with A'^tot. > -^obs. it is not possible to make firm predictions 
for the observable signatures in our horizon patch but one is led to consider probabilities 
of different signatures. It turns out that especially in non-Gaussian models the differences 
between the entire inflating patch, and of a random patch the size of our observable patch, 
can be significant. This is the case even if the horizon exit of the largest observable modes 
would be preceded by just a few e-foldings of inflation. 

Since the fundamental model of inflation should be confronted with theories of high- 
energy physics, it is of cruciai importance to develop a solid understanding of its relation to 
the effective description of the last A^obs. e-foldings, which in turn can be directly confronted 
with the observations made in our horizon patch. In this work we will make systematic 
progress towards this goal. We will establish novel relations between the observable signa- 
tures and the form of the underlying inflationary physics. Our model independent approach 
opens up interesting new possibilities to address the naturalness of different non-Gaussian 
signatures. By making use of these findings, we will establish new connections between 
large and small non-Gaussianity and the structure of inflationary physics. 

2. The curvature perturbation 

We concentrate on the class of generic single source models where the dominant contribu- 
tion to curvature perturbation is due to a single scalar field, known as single-source inflation 
||5|. This field can be an isocurvature degree of freedom during infiation, e.g. the curvaton 
scenario |^ or modulated reheating ^. Such models can easily generate observable 
non-Gaussianity. In the presence of isocurvature directions, the long-wavelength fluctua- 
tions can not be removed by merely shifting the time coordinate like in the adiabatic single 
field case Consequently, the long-wavelength modes can have interesting observational 
consequences. 
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We assume the non-Gaussianity is generateci on superhorizon scales and neglect ali 
contributions from subhorizon modes. While some non-Gaussianity is in general also pro- 
duced by subhorizon physics, these contributions are slow roll suppressed in models with 
canonical kinetic terms. In this class of models the leading contribution to observable 
non-Gaussianities is therefore generated by superhorizon physics. 

Neglecting the subleading subhorizon contributions, the curvature perturbation ^ can 
be expressed in the locai form 

c(x) = cg(x) + ^/Sl (cè(x) - m + (^)%&LCà(x) + [l)'h'^L (cè(x) - m + ■■■ 

(2.1) 

Here Cg is a Gaussian field and the tree-level coefficients f^^, g^^ etc. are constants. We 
are neglecting a possible scale dependence of the non-linearity parameters, as discussed for 
example in [|lO|. 

The locai Ansatz amounts to assuming that C is an analytic function of a single Gaus- 



sian inhomogeneous field. To make this manifest, we can equivalently rewrite (2.1) as 



n=l 



where ero denotes a classical homogeneous solution of equations of motion (and we under- 
stand that we remove from the even powers in the previous formula the averages over the 
entire space, to make the average of C over ali space vanishing). This corresponds to the 



familiar 6N expression |11| if 5cr(x) = (Q{x)/N'{ao) are chosen to represent fluctuations 
on a spatially fiat hypersurface soon after the horizon exit. N{a) then measures the num- 
ber of e-foldings from that hypersurface to a uniform energy hypersurface at which C is 
frozen to a Constant value (which must exist provided that ali isocurvature perturbations 
decay) . 

To set the notation: the spectrum, bispectrum and trispectrum are parameterized by 

(CkiCk^) = {27rf6{Y^ki)Po{h) , (2.3) 

(CkiCk.Cka) = {27rfS{Y,^i){f^Mh)Po{k2)+peTms.) , 

(CkiCk,Ck3Ck4) = {27rf6iY^ki){g^^Mki)Po{k2)Po{k3) + 

+ T^Mki)Poik2)Poik3)Poi\ki + k2\) + peTms.) . 

Since we are considering single source scenarios, the tree-level parameters /^l ^^'^ ''"nl 
related by r^^ = (6/5)2 (/«J^ ^ |1. 



2.1 Long w^avelength modes 

Let US consider the impact of long-wavelength fluctuations of the field a in the expression for 
the curvature perturbation (p.2|). We denote by ero the classical homogeneous solution which 
corresponds to the spatial average of o"(x) over the entire inflating patch. Fluctuations 
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around the global background can be divided into long and short wavelength components 
with respect to the horizon scale of the observable universe fcobs. = «obs.^^obs. 

S^i^) = [ 7^ e'^--6a{q) + j ^ é^--5a{ci) ^ 5a,{^) + 5aU^) . (2.4) 

To first order in perturbations the long- and short-wavelength modes are uncorrelated, 
(5crL(k)(5(Ts(k')) = 0, and they both have vanishing ensemble averages, ((5o"l) = {Sas) = 0. 
The ensemble averages are assumed to coincide with the spatial averages computed over 
the entire inflating patch, hence (it(x)) = (Tq. 

Spatial averages computed over regions smaller than the inflating patch in general 
diff'er from the ensemble averages. This is due to the long-wavelength modes. For example, 
the average of óaj^ over a spherical patch with radius k~^^ and origin at xq reads 

((5CTL(x))obs.,xo = TT— I <^^^(.Khs ~ |x-xo|)5aL(x) (2.5) 

l^obs. J 



- <5o-l(xo) . 

In order to pass from second to fast line, we have assumed a scale invariant form for the 
mode functions, 5a{q) oc g"^/^. 

In a similar way one finds (5o'£(x))obs.,xo — ^^^lÌ^o), so that the curvature perturbation 
( |2.1| ) in a patch the size of our observable universe k~^^ can be written as 



n=l 



ni 



Here 5cjl(xo) is a Constant in (...)obs.,xo (but not in (...) since it depends on the position 
Xq). We have adopted a notation where ^obs. denotes a quantity X evaluated in a patch 
of the size of our observable universe. 

The results depend on the location xq of the patch as the contribution of long- 
wavelength modes is difFerent in different regions of the inflating space. This position 
dependence limits our ability to make precise predictions for observables in our horizon 
patch 10 if A^tot. > -^obs.- Instead of assuming that we would occupy a "typical" loca- 



tion as in |14] we follow a top down approach and systematically relate the probabilities 
for difFerent observational configurations in our patch to the form of the underlying fun- 
damental inflationary physics. Given an inflationary model with initial conditions set at 
the beginning of inflation, we can make precise predictions only for spatial averages over 
the entire inflating patch. These can signiflcantly differ from the averages in our horizon 
patch due to the unknown long-wavelength contributions. However, given a model we can 
unambiguously compute the probabilities for different observational conflgurations in our 
horizon patch and thus work out the possible observational imprints of the model. Going in 
the opposite direction, we can also convert the observational Information into constraints 
on the fundamental high-energy physics controlling the evolution of the entire inflating 
patch. 
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In this work we will systematically explore the consequences of the long-wavelength 
modulations. We will establish novel connections between the effective inflationary physics 
describing our observable patch, and the fundamental inflationary physics which dictates 
the statistics of perturbations over ali of the inflating space and directly reflects the un- 
der lying high-energy physics. 



3. Probability distributions of locai observables 

If inflation lasted longer than A'obs. ~ 60 e-foldings, our observable universe constitutes 
an exponentially small fraction of the entire inflating space. Patches of our horizon size 
(aobs.-f^obs.)"^ located in different places of the universe can then have different physical 
properties. This is a consequence of the long-wavelength fluctuations which generate a 
random modulation in the locai expansion history. The modulation is proportional to the 
number of e-foldings before the horizon exit of our observable region 

iVin. ^ In , (3.1) 

where to stands for the beginning of inflation. The longer the duration of inflation, the 
more the unknown long-wavelength contributions limit our ability to directly probe the 
underlying high energy physics reflected in the statistics of inflationary perturbations over 
the entire inflating space. 

3.1 Spectrum, bispectrum and trispectrum 

We assume that the fluctuations at horizon crossing 6a = Cg/^' in ( p.2| ) are Gaussian so 
that ali correlators can be expressed in terms of the two-point function 

{Sa{k)6a(k')) = (2^)35(k + k')Ps^{k) = {27Tf6{k + k')-^n<x • (3.2) 
For the canonical models of inflation we are concentrating on, 

■Psa ^ {H/2-Kf and the 

dependence of H{t, a) on the fleld a is slow roll suppressed. In the following we will 
neglect this dependence so that the two point function of 6a computed over an arbitrary 
subregion of the inflating patch is given by the ensemble average (|3.2| ) . In what follows we 
will also ignore the slow roll suppressed scale dependence of Vsa- 



Using (2^), the spectrum of curvature perturbations in a patch the size of our observ- 



able universe is given by 

Vobs.=Vo(^l + ^f^^N'{ao)óaU^o)^ , (3.3) 

to first order in the long-wavelength field. Vq = N''^{fTo)Vsa and = 5A^"((To)/6iV'2(cro) 
denote the ensemble expectation values of the spectrum and the bispectrum amplitude, 
which coincide with spatial averages over the entire inflating patch. 

Hence, the spectrum "Pobs. of fluctuations in a horizon size patch depends on the lo- 
cation xq of the patch through the Gaussian field (5(Tl(xo). This gives rise to a landscape 
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picture where the amplitude of curvature perturbations measurable locally on patches of 
size (oobs.-ffobs.)"^ fluctuates around the global average Vq according to a Gaussian prob- 
abihty distribution 



PiVoi^s.) = (2vra^)-^/^ exp ( ) . (3.4) 



,-,.,,.^2.^-1/2 

Here the variance is given by 



Mn. . (3.5) 



In a similar way, using equation ( |2.6| ) and working to first order in Jcjl , we can calculate 
the probabiUty distributions for f^', and g^^' measuring the non-Gaussian statistica 
in a horizon size patch. The results are given by 

expf-^^fc^^^ 



27rcrf 

/NL 

exp ( — — j 

^'tó') = ^ , "^^"^ ^ , = ( T^NL - ) T'oA^.n. , (3.7) 

where /i^^ defined as the fourth order coefficient in the expansion of the curvature 
perturbation (p.l|). The non-hnearity parameter in each patch is given by = 
(6/5)^ (/§L^)^ at scales corresponding to the patch size. This is a result of the single source 
forni of the curvature perturbation Cobs.( P-6D - smaher scales the relation may become 
violated by the loop corrections as discussed in [p^ ]. 

Notice that the dependence of inflationary observables on the location of the patch is 
closely associated to semiclassical loop contributions to observable quantities, as we inves- 
tigated in ||2|, 15] in collaboration with D. Wands. The difference between the curvature 



perturbation Co over the entire inflating patch and the curvature perturbation ^obs. in our 
observable patch is due to the long-wavelength modes aoHo < k < aobs.-^obs. generated 
before our horizon exit. Integrating over the unobservable long wavelength modes to find 
Cobs. ( |2.6D amounts to computing radiative corrections to Co at the scale /cobs. = Oobs.-f^obs.- 
As explained in the radiative corrections to an n-point function of curvature pertur- 

bation consist of both the usuai loop corrections, where internai momenta are integrated 
over, and of soft limits of higher order n-point functions where external momenta become 
unobservably small. The soft contributions generate variances for the n-point functions of 
Cobs. causing the properties of a patch the size of our observable horizon to depend on its 
location within the entire inflating space. It is precisely the ramifications of this effect that 
we are investigating in the current work. 

As also explained in the next to leading order long-wavelength corrections to 

"Pobs. ) /nl^' ) T'Nif ' ) 5nl ' associated with higher powers of Vq are subdominant with respect 
to the leading order corrections as long as the bare non-linearity parameters f^^, g^^, 
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etc. are not extremely large. In this work we are not considering such extreme values 
and our first order analysis of the long-wavelength corrections is therefore justified and self 
consistent. 



3.2 How large are the variances? 

The variances generated by the long-wavelength fluctuations are uninterestingly small if 
there was basically no inflation before the horizon exit of the largest observable modes, 
A^in. < and if the ensemble expectation values of the non-Gaussian amplitudes are 

small, I/nlI ^ ^(1)) I^nlI ~ (/nl)^- "^^^ more one deviates from either of these conditions, 
the bigger the variances become. 

The growing variances affect the comparison of inflationary models and cosmological 
observations in two different ways. As the variances of the theoretical predictions for the 
observables get bigger than the observational sensitivity, the bounds on model parameters 
get weakened, and projection of observational constraints on the space of model parameters 
may also become non-trivial. Comparing the variance cip of the spectrum Vohs. ) given by 
equation ( |3.5D to the 1-a error AV/V — 0.1 in WMAP 7-years parameter fits, we find 
that the theoretical uncertainty dominates over the observational inaccuracy if 

aT,>AV ^ \lm>(!^] . (3.8) 



50 ~ V300, 

Here we have set = 2.44 x 10^^. For inflationary modes featuring a few hundreds of e- 
foldings, iVin. > O(102), the variance becomes comparable to the observational accuracy for 
I/nl! ~ O{10). This has the effect of broadening the region of parameter space compatible 
with observations. The growing variance increases the class of models with different super- 
horizon scale properties, and underlying inflationary physics, but degenerate predictions 
on observable scales. 

The variance of the spectrum can dominate over the expectation value, cr-p > Vq, if 
the inflationary model features extreme non-Gaussianities, e.g. I/^lI ^ lO^i and/or a very 
long period of inflation A'^in. ^ 10^. In both cases one should go beyond first order in the 
long-wavelength perturbations (5ctl. A more subtle issue is that the validity of the entire 
semiclassical approach, justifying the use of the expansion (|2.2[) , may become questionable 
in the regime where fluctuations grow large: this is the regime investigated for example in 
Q in the case of the curvaton scenario. In this work, we restrict ourselves to the regime 
of small perturbations, tractable by a first order analysis in (5ctl, where the semiclassical 
approach should provide an accurate description of the system. Consequently, we will not 
encounter situations where the variance a-p would dominate over the global background 

For the bispectrum the situation is different. Comparing the variances of the spectrum 
(|3.4|) and bispectrum (^), we obtain the relation 
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Clearly, for models with |(7nlI ^ (/nl)^ ^he relative variances are of equal magnitude. Using 
the expected 1-a sensitivity of Planck, A/nl = 5 |16], we then find that the theoretical 
uncertainty for f^- becomes comparable or dominates over the observational accuracy, 
f^/NL ^ ^/nl, for 

I/SlI>5(^)"' . (3.10) 

As we will discuss below, this can happen in the regime a^p < AV where the theoretical 
variance of the spectrum is stili small compared to the observational accuracy, and we may 
set T^obs. — "Po to a reasonable accuracy. 

On the other hand, for models with (/nl)^' ^he variance of the bispectrum can 

become very large even if our horizon exit was preceded by just a few e-foldings of inflation. 
It may even dominate over the background value, cj/^j^ > I/nlI' effectively screening the 
underlying inflationary physics which determines the global structure of the inflating space 
and reflects the fundamental high energy physics behind inflation. While the global average 
/^L represents the most probable value also in patches the size of our observable universe, 
the probability for observing a drastically different bispectrum amplitude becomes 
significant as the variance CFf^^^ grows large. This is illustrated in Fig. ^ for the specific 
example where /^l = globally but the model may stili result a detectable bispectrum 
in our observable patch. 




Figure 1: The probability to find 10 < f^l^- 80 in our observable part of the universe for 
different values of assuming that the bispectrum averaged over the entire inflating patch 

vanishes /^^ = 0. The probability is plotted against, N-m., the number of e-foldings from the 
beginning of inflation up to the horizon exit of our observable patch. The three curves correspond 
to g^L — (red, rightmost), g^L = 5 x 10'' (green, middle) and g^L — (blue, leftmost), and 
we have set Vo = 2.44 x 10~^. 

The variance of the trispectrum may also be large. As we are considering single 
source models, the observable calculated using ( |2.6|) is, at the scale /cobs. = Oobs.-ffobs. 
corresponding to the tree-level in the observable universe, completely determined by /§jf' 
through the relation = (6/5)^ (Z^^*^')^. Its variance is then given by 

= 4 = Hr^L?^, ■ (3T1) 
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It dominates over the background, cr'^^ > (t^l)^, under the same conditions as cr'j^-^ ^ 
(/nl)^- Looking at the equation ( p.TD , we also find that the variance of g'NL niay become 
large, 0"^^^^ > (S'nl)^' if the amphtude of the 5-pt function is large, I/i^lI I^nlI- However, 
we wiU not consider this last possibihty in this work. 



4. Implications on inflationary physics 

The cumulation of long-wavelength fluctuations can generate significant variances in the 
smah scale properties of primordial perturbations across the inflating space. If inflation 
lasted much longer than the A'obs. ~ 60 e-foldings corresponding to the horizon exit of 
the largest observable modes, this leads to a landscape picture where the properties of 
our observable patch depend on its location within the entire inflating space. This is 
qualitatively similar to the special case of the curvaton web discussed in 

The superhorizon variation is not measurable and the observable perturbations can be 
described by an effective model of inflation, covering only the last A'^obs. e-foldings. The 
effective model may however significantly differ from the fundamental model for the entire 
inflating epoch. The Lagrangian of the fundamental model is directly determined by the 
under lying high-energy physics which may also provide a naturai range of initial conditions. 
In the effective model, the field values in the Lagrangian are shifted by random long- 
wavelength fluctuations. This may significantly affect the apparent form and predictions 
of the theory in the presence of isocurvature directions. 

The ensemble averages calculable from the fundamental model correspond to the most 
likely predictions also in patches smaller than the inflating volume, at least in the per- 
turbative regime we are considering here. The variances however make it possible to find 
substantial differences between the locai and global properties of primordial perturbations. 
It is therefore of key interest to carefully explore the relation between the observable locai 
quantities and the unobservable global inflationary perturbations which reflect the under- 
lying high energy physics. In particular, the differences in non-Gaussian statistics over the 
entire inflating patch and our observable universe may easily become large with interesting 
consequences for observable quantities. 

4.1 Small primordial non-Gaussianity 

If Planck would not detect primordial non-Gaussianity, \f^\ ^ 10, then single-field slow- 
roU inflation would be favoured as the minimal scenario consistent with observations. How- 
ever, nature might not have chosen the minimal setup. A small bispectrum |/§jf| ^ 10 
could also arise from models with signiflcant non-Gaussianities in the form of trispectrum 
I^nl'I > ifm!' f (see 0, |8| for examples of such scenarios). The signal would remain 



unobserved by Planck if \g^^' \ < 10^ but could be revealed by future large scale surveys. 

Using the landscape picture we introduced above, we can obtain novel insights on the 
naturalness of such models. The variance ( |3.6D of the locally observable bispectrum am- 
phtude grows as cr/^j^ oc ^^L-^r/^ assuming that (/^l)^ ^ 5nl- Therefore, even if 
the global bispectrum would be small I/^lI < 10) probability for obtaining an unob- 
servably small bispectrum amplitude \f^'\ ^ 10 in our patch decreases as grows or 
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the inflationary epoch becomes longer. This is illustrateci in the Fig. |2| which depicts the 
probability for \f^^j^-\ < 10 as a function of /^l, g^i^ and A^in., characterizing the fundamen- 
tal inflationary model. For the parameter range shown in the figure, the variance of the 
spectrum is negligible compared to the observational accuracy a-p <^ AV and we have thus 
set Vq ~ "Pobs. = 2.44 X 10~^. The variance of trispectrum is also small so that g^^' ~ g^j^ 
to a good precision. 



^Nl Nin.= 1 gNL Ni„.= 10 gNL Ni„.= 100 




Figure 2: Probability to obtain a non-detectable bispectrum I/^l*^ ] < 10 as a function of the 
ensemble expectation values /^l .9nl characterizing the underlying fundamental inflationary 
model. Three values of Mn. are shown corrcsponding to different total durations of the inflationary 
epoch. 

The figure shows that models with I/^lI ^ 10 and I^'^lI ^ 10^ have a high probability 
P > 0.9 of producing unobservably small bispectrum in our horizon patch. The prob- 
ability is not sensitive to the duration of the inflationary epoch. We thus conclude (within 
the limits of validity of our perturbative treatment) that a small bispectrum I/nl' I < 10 
is a naturai outcome in inflationary models with a Lagrangian leading to I/^lI ^ 10 

b^Ll < 10^- 

On the other hand, Fig. Q also shows that as the trispectrum grows observable I^^lI — 
I^nl'I ^ 10"^ the probability for an unobservable bispectrum I/nl' I < 10 becomes smaller 
but not negligible. It is also interesting to note that an unobservable bispectrum can 
be generated even in models with a large global I/^lI ^ 10 even though the probability 
for this outcome is not higher than a few tens per cent at most. In the regime of large 
trispectrum j^^Ll ^ IO'' we also note that the probabilities become sensitively dependent 
on the amount of inflation A^in. before our horizon exit. For example, the Fig. |2| shows 
that the probability to find \f^'\ < 10 and I^nl*'! ~ 10^ becomes very small for A^'in, » 10 
while it may be substantially larger for A^in. < 10. Therefore, to generate a signature 
I/nl^ I ^ 10; I^nl I ^ 10^, inflation should not have lasted much longer than the observable 
-Aobs. e-foldings which represents a non-trivial tuning of model parameters. 

We emphasize that our approach is model independent, apart from our restriction to 



the locai form of non-Gaussianity (2T). We therefore cannot address the possible fine- 
tuning associated to realizing specific values for /^^ ^^à. g^^ in a concrete inflationary 
setup. In contrast, our approach reveals that, independently of the details of inflationary 
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physics, certain configurations /nljS'nl general fail to translate into similar observa- 
tional configurations /^if ' , S'nl ' unless non-trivial constraints are placed on the duration of 
inflationary epodi. We will address this model-independent tuning in more detail below. 

4.2 Large primordial non-Gaussianity 

In the regime of detectable bispectrum by Planck, 1/^1^1 ^ variance due to long 

wavelength modes can have interesting effects. As the variance grows bigger than the 1-a 
accuracy of Planck, (T/j^^ ^ ^/nl = 5, the projection of observational constraints on to 
the model parameters /^l ^^à. g^i^ becomes non-trivial. 

This is illustrated in Fig. ^. The figure depicts the relation between the observable f^, 
and the theoretical parameters /^l ^^'^ 9nl directly connected to details of the underlying 
inflationary physics. The variance (Jf^^ grows along with the number of e-foldings A^in. 
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Figure 3: The observable bispectrum /^l ' plotted against the theoretical expectation value /^^ 
which reflects the underlying inflationary physics. The shaded area corresponds to the observational 
1-a uncertainty of Planck A/nl = 5. The curves in the upper panel bound the theoretical 2- a 
regions f^^- = /^l ^ '^'^ìnl the bispectrum amplitude f^^- in our observable patch. The 
lower panel shows the corresponding 1-a regions /^if' — /nl ^ '^/nl- The three curves depict 
9nl = (black, middle), = (/Sl)^ (^'^d, innermost) and = 10(/^l)^ (blue, outermost). 
The results are shown for three choices of Ni^,, the number of e-foldings before the horizon exit 
of largest observable scales. The variance of the spectrum is smaller than the observational error 
c-p/'Po < 0.1 in the entire regime shown in the plots. 

preceding the horizon exit of largest observable modes. For I^^lI ~ (/nl)^' " region 
of can grow so large that a detection of would not suffice to place tight constraints 



- 11 - 



on /^L- Consequently, the observable value /^if' may considerably deviate from the global 
average f^^. The same would also be true for the trispectrum g^^' if I/i^lI ^ I^nlI' 
the 5-pt function would be large. Here we however do not consider this possibiUty further 
but restrict ourselves to models where \h^^\ < \g^i^\, and consequently g^£' ~ g^^^. 

We recap that /^^ corresponds to the tree-level amplitude of the bispectrum computed 
for the inflationary model describing the evolution of the entire inflating patch. The forni 
of this fundamental inflationary model is eventually dictated by the underlying high-energy 
physics and we denote its Lagrangian schematically by Cq. On the other hand, the prop- 
erties of our observable patch can be parameterized by an effective model £obs. describing 
only the last A'obs. ~ 60 e-foldings of inflation in our patch. Observables such as then 
denote quantities calculable from the effective model >Cobs.- In the semiclassical approach 
that we are utilizing here, the field values of the effective Lagrangian are shifted by the 
random long-wavelength contributions, as is clear from equation ( p.6D . From the point of 
view of the effective model the shifts correspond to random changes of initial conditions 
if isocurvature directions are present during inflation. These can significantly affect the 
predictions, and especially the non-Gaussian statistics, generating substantial differences 
in the apparent form of the fundamental model Cq and the effective description of our 
patch £obs.- This has interesting ramifications as only the fundamental model Co is the 
one directly reflecting the high energy physics behind the inflationary epoch whereas the 
effective model >Cobs. is the one directly constrained by observations. 

As the Fig. ^ illustrates, the variance cr/^^ can be large compared to the observational 
1-fj error A/nl not only in models with I^^lI > (/nl)^ ^ut also for I^^lI ^ (/nl)^ /nl 
is large enough. As a curiosity, it is interesting to note that the variance is accidentally 
suppressed around g^i^ ~ (/nl)^ "^^^ to cancellation^ of the two terms in (|3.6D . In this 
regime the predictions of the model are not significantly affected by the duration of inflation 
but the theoretical signature g^^^ ~ (/nl)^ i^ generically reflected with a high probability 
by the observables g^^' and f^'- This behaviour is illustrated in Fig. |^. 

When the variance a/^j^ grows comparable to I/^lI the landscape effects become even 
more pronounced. In this regime the observational signature becomes dominated by the 
long-wavelength modes which may effectively screen the originai signal /^lj reflecting the 
structure of the underlying inflationary physics. The locally measurable /^jf' then becomes 
strongly dependent on the location within the inflating patch and /§jf' may significantly 
deviate from the ensemble expectation value Z^^. This behaviour can be observed in Fig. ^ 
which depicts the probability of generating an observable bispectrum as a function of the 
model parameters /^l ^nd g^i^. Indeed, in the regime I^nlI ^ (/nl)^ where fj^^^ ^ (/nl)^ 
the probability for obtaining 10 < < 80 remains non-negligible even for /^^ ^ 100- 
The boundaries of the 60% probability contours in the /nl" direction however do not 
significantly differ from the bounds 10 < /^if ' < 80 unless the inflationary period becomes 
very long, Ni^i > 10^. The slow tightening of the upper bounds on /^^ ^-s a function of A'^in., 
as observed in the figure, is due to the growth of the variance Cf^^^ ( |3.6D . In the regime 

^The vanishing of af^-^ for gr^L = 4/3(/Sl)^ (the proportionality factor depends on the Ansatz for the 
curvature perturbation, in our case eq. (2.1)) is an artefact of neglecting the slow roU corrections; in general 
they would generate a small (Jf-^^ even at — 4/3(/Sl)^- 
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Figure 4: The countours depict the probabihty to obtam a bispectrum in the range 55 < /^if ' < 65 
as a function of the model parameters /^Jl f^nd g^L- The resuhs are shown for three choices of iVin., 
the number of e-foldings before the horizon exit of largest observable scales. The variance ct/^l 
is accidentally suppressed for g^^^^ ~ (/nl)^- This enhances the probabihty to produce larga 
close to the ensemble average under lying inflationary physics. 




Figure 5: Probabihty to generate a bispectrum in the observable range 10 < |/Nif | < 80 as a func- 
tion of the ensemble expectation values /^l ^^.d Three values of N-m. are shown corresponding 
to different total durations of the inflationary epoch. 

5nl ~ (/nl)^ variance is suppressed which explains the different behaviour observed 
for these values in Fig. |^. 

For even longer periods of inflation Ni^^^ ^ 10^ the effects are expected to become even 
more significant. The analysis of this regime is however beyond the scope of the current 
work as our perturbative approach can not be directly applied to this regime without 
carefully addressing convergence issues. 

4.3 Non-Gaussian signatures and fine-tuning 

If the variances grow large, the observational configurations become sensitively dependent 
to A'^in., the amount of inflation before our horizon exit. For the bispectrum this happens 
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when I^nlI ^ (/nl)^ '^^^ clearly observed in Figs. ^ and |5[ 

Indeed, our findings show that observable patterns {/nl'',5nl'} with IgN^L'l ^ (/nl' )^ 
directly reflect the under lying theoretical signature {/nl'5nl} o^^y ^hat signature was 
reahzed just before the horizon exit of the observable scales. Consequently, as {/nl'^nlÌ 
characterise correlators over the entire inflating space, this impHes that inflation should 
not have lasted much longer than the observable A^^^bs. e-foldings. If inflation lasted much 
longer, A^'in. ^ 1, the rapidly growing variance of the bispectrum 0"/^^ effectively screens 
the originai signature {/nl'S'nlI- "^^^ probability of deviating from /^^ then be- 
comes significant and the observational configuration {/Nif',5'NL } longer reflects the 
theoretical signature {/nl'^nl} dictated by the underlying inflationary physics. Related 



considerations were also made in 1171. 



There is no a priori reason why the inflationary epoch should not have lasted mudi 
longer than A/'obs. e-foldings, corresponding to the largest observable modes in the current 
universe. Consequently, requiring A'in. to be small is a non-trivial constraint indicating 
generic tuning associated to configurations with farge g^^' ■ The range of e-foldings A'^in. 
after which an originai signature Is'nlI ^ (/nl)^ S^*^ screened by the long-wavelength 
modes is depicted in Fig. |^ as a function of g^i^. Note that since we assume the amplitude 
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Figure 6: Configurations I^nlI ^ (/nl)^ appear tuned as they become strongly modulated by 
long-wavelength modes if the horizon exit of the largest observable scales was preceded by even a 
few e-foldings of inflation. The dashed line depicts the vaine chosen for /^^ i the amplitude of the 
bispectrum over the entire inflating patch which reflects the underlying inflationary physics. 



of the 5-pt function is not exceptionally large, 

,2 



1/7° I 



< 



the variance of g^^' is 



negligible, cr'g^^ <^ (Snl)^' range of parameters we are considering. Therefore to a 



5nl- 



good precision we can approximate g^L — 

Fig. ^ shows that for signatures with ^f^^^- 10^ and (f^')"^ ^ S'nL ' bispecturm 
is dominated by the long-wavelength contributions already if A^in. ^1- To translate a 
theoretical configuration {/nl'S'nl} with I^^lI ~ ^° corresponding observational 
signature {/nl^ ^^nl }' parameters ^-i^d g^^ should thus be tuned to the desired 
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values less than one e-folding bcforc the horizon exit of the obscrvablc mode. If the values 
/^L ^^rs set at an earher stage, the observable has a high probabihty to 

significantly differ from the theoretical expectation Z^^- This suggests that such config- 
urations are not very naturai - unless the model under consideration naturally generates 
-^tot. ~ -^obs. e-foldings of inflation. For \g^' \ > 10^ the corresponding range of e-foldings 
is given by A^in. < 10 indicating a somewhat milder tuning. For I^^l ' I ^ the constraint 
is relaxcd to A'^in. < 100. 

We stress again that our discussion of the fine-tuning addresses only the probabilities 
of finding observational confìgurations {/Nifs^Nx'} in our horizon patch given that our 
horizon exit was preceded by some e-foldings of inflation. The naturalness of realising 
certain values of /nlj^nl specifìc inflationary models is another, model dependent issue 
which may give rise to further tuning, or favour specific signatures. Our general approach 
instead addresses the way how the theoretical confìgurations g^^^ translate into 

measurable properties in our observable patch. As we have shown, for the class of models 
with \g^i^\ » (/nl)^ these considerations suggest that converting the vaine /^l i^^o the 
corresponding observational value ~ /^^ becomes unlikely unless the horizon exit of 
our observable patch occurs soon after the beginning of inflation. Imposing this requirement 
indicates a generic source of non-trivial tuning in models with I^^lI ^ (/nl)^- 

5. Extension and applications of our results 

5.1 Variance of the power spectrum in multi-source scenarios 

In this paper we have focused on single-source scenarios, and shown under which circum- 
stances there may be a significant variance of observable quantities when making observa- 
tions in a patch smaller than the total inflated volume. Redoing this analysis in multi-source 
scenarios goes beyond the scope of this work, but our formalism makes it possible to do 
so, and here we give one small example, thereby showing that our considerations remain 
interesting in multi-source cases. 

In a multi-source scenario, we must take into account the effect of ali relevant long 
wavelength fluctuations on to the short wavelength modcs. For the power spectrum, the 
lowest order correction comes from the influence of the linear long wavelength modes on 
to the linear short wavelength modes 



where summation convention is used for the field indices, A,B,---, over ali fìeld compo- 
nents. 

At tree level, and assuming that each component of the field perturbations are uncor- 
related, we may use 



Cobs. = KScT^ + N^sScT^SaE + ■■■ 



(5.1) 



a ì 




(5.2) 
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Hence 



(Cobs.,k,Cobs.,k.) = A^^iV^(Kki'^<k.) + (A^°^Bc{'^<k/<k.)<5<^L + 1 Perm.) (5.3) 
^obs.=7^o(^l + 2^^Li.j , (5.4) 
= ^rl^VlN,^. . (5.5) 
Applying the Suyama-Yamaguchi (SY) inequality r^^ ^ EU' follows that 

Hence for a given value of /^l, the variance of the power spectrum is larger than in single- 
source scenarios, showing that our consideration about the difference between observables 
in the entire inflated patch and our observable patch may become even more important in 
multi-source scenarios. 

5.2 What will we observe if the single-source equality is broken by loops? 

Since we have assumed in most of this paper that we are in a single-source scenario, it 
is interesting to consider how we can observationally check whether this is the case. The 
weh known test is the equahty T^jf' = (6/§l'^'/5)^ which however is only vahd at tree-level 
@, |0|. In what follows we are mainly considering quantities in our observable patch and 
for brevity we omit the labels where no risk for confusion arises and simply denote 
/nl = Z^if' etc. Loop effects can most easily break the equality when (^nl is large, a 
scenario we have extensively studied in this paper, making this particularly relevant to 
study here. This would cause a deviation from the equality when considering loop corrected 
quantities at scales smaller than our horizon patch. However there is a generalised equality 
between radiatively corrected quantities which include both loop corrections where internai 
momenta are integrated over and soft corrections where external momenta become small. 
This new equality is valid to ali orders in radiative corrections and it provides a new test 



of inflation, as shown in jl5|. 



A priori, in order to see a breaking of the equality, tnl is the only parameter which 
we need to observe, regardless of whether it is broken by loops in a single-source scenario, 
or by multi-source effects. Here we will show that in practice if it is broken by loops, 
we are essentially guaranteed to see ali three non-linearity parameters. The calculation is 
restricted to the one loop level. 

To be concrete, we assume that Planck can detect non-linearity parameters with 2 — a 



confidence if /nl > 10,rNL > 10'', ^nl > 10 |^ 



Assuming that the SY inequality is always valid |23, 2C, 24, 25 j, to see a breaking of the 
equality with Planck we have to observe 10'^ < tnl ^ 10"^ where the upper bound comes 
from the current observational constraint. In this case 30 < |/nl| ^ 10^ if the equality 
holds and we would clearly observe /nl as well. The key observation is that in order for 
/nl not to be seen, |/nl| < 10, we need ^nl = tnl/(6/nl/5)^ > lO^/lO^ ~ 10, so the 
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equality needs to be strongly broken in this case. From [15| we find at one-loop level 

A„..(^)'te) .,+4(f!t)' V, (5.7) 

\"/ \ /nL / trec+loop \JNL/tree 

where the second term on the right hand side of the last equahty is responsible for breaking 
the equahty. In order to get ^nl ^ 10 with |/nl| ^ 10 we should then have |(7nl| ^ 
3 X 10^ ^ tree-level and hence the "normal" hierarchy between (Jnl and /nl must 

also be strongly broken. A breaking of the equahty must therefore come through an 
observably large ^nl, which is dose to the current upper hmit. Such a large tree-level ^nl 
does not create problems with a loop correction to either /nl or the power spectrum, and 
the perturbative expansion in loops remains valid [p^ j. 

Finally given that we require such a large value of what is the probability that 

is too small to be seen in our Hubble volume value? The answer can be seen in Fig. |2|, 
for I^NL^'I I^nlI ^ 10^ probability of obtaining I/nl^'I < 10 is suppressed below about 
50% already for A^in. 1 and it rapidly decreases as A^in. grows. We may thus conclude 
that a breaking of the SY equality in single-source inflation is likely to be accompanied 
by a measurement of ali three non-linearity parameters. It should then also be possible 
to observe our equality rjS|^' = (6/5) ^(Z^'^')^ between the radiatively corrected quantities 

and /nl^' which contain not only the loop corrections but also the corrections from 



soft modes P, 15]. 



6. Conclusions 

Primordial perturbations with wavelengths greater than the observable universe contribute 
to background quantities in our observable patch. The long-wavelength contributions vary 
across the inflating patch and amount to random shifts of the locai field expectation values. 
This leads to a landscape picture where the properties of our observable patch depend on 
its location. The observable primordial perturbations constrain the form of the effective 
model for the last A^^bs. e-foldings of inflation in our observable patch. Being affected by 
random long-wavelength contributions the effective model may considerably differ from the 
fundamental model of inflation which describes the entire inflating patch and whose form 
is dictated by the underlying high energy physics. 

In this work we have shown that if multiple fields are present during inflation, the long- 
wavelength contributions can be significant even if our horizon exit would be preceded by 
only a few e-foldings of inflation. Long-wavelength modes generate shifts along isocurvature 
directions which correspond to non-trivial changes of initial conditions in the effective model 
and can greatly affect the resulting perturbations, especially the non-Gaussian statistics. In 
particular, the impact of the long-wavelength modes can not be deduced from the effective 
model itself but they reflect the form of the underlying fundamental model of inflation. This 
is in sharp contrast to pure single field slow roU models where the long-wavelength effects 
are slow roll suppressed and amount to locai time shifts expressible in terms of spectral 
indices. As the simplest example of multiple field models, we have in this work concentrated 
on long-wavelength effects in the generic class of single source models with the locai type 
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of non-Gaussianity. Many theoretically interesting models, for example the curvaton and 
modulated reheating scenarios in the hmit of neghgible inflaton perturbations, fall into this 
class of models where perturbations are generated by a single inhomogeneous field, despite 
the presence of several light fields during inflation. 

We have found novel general relations between the primordial perturbations in our 
observable part of the universe and the inflationary physics controlling the entire inflating 
patch. We have investigated the consequences of both a detection and non-detection of 
the primordial bispectrum /^if' by Planck. Our results demonstrate that an eventual non- 
detection would not strictly rule out models with a large bispectrum. Indeed if inflation 
lasted longer than the observable A'obs. e-foldings, models with a bispectrum up to /^^ ~ ^0 
when averaged over the entire inflating patch, can have a probability up to few tens of per 
cents for yielding l/^if l ^ 10 in our observable patch. However, this is not a very generic 
feature but requires a large trispectrum amplitude I^^lI ^ ^ vaine of falling 

to a specific narrow range. 

In case /nl is detected, the long-wavelength effects quite generically have interesting 
implications for the interpretation of the observation. If the inflationary model predicts 
5nl ^ (/nl)^ across the entire inflating patch, we find that the observable bispectrum Z^^^- 
has a high probability to deviate from the global mean /^^ even if there were only a few 
tens of e-foldings of inflation before our horizon exit. As we have shown in this work, the 
observational constraints on therefore translate in a highly non-trivial way on the 
parameters f^^, g^^, characterizing the under lying inflationary physics. 

We have also shown that inflationary models predicting a hierarchy I^'^lI ^ (/nl)^ 
imply an inherent tuning of the total duration of inflation irrespective of details of the 
inflationary physics. This is because the variance of the bispectrum easily becomes very 
large after a short period of inflation making it unlikely to find an observable vaine /§if ' 
which would reflect the global mean f^^. For I^^lI ~ ^^^^ happens if our horizon 
exit was preceded by just one e-folding of inflation and for I^^lI ~ after about ten 
e-foldings. To translate a theoretical configuration Is'nlI ^ (/nl)^ ^° corresponding 
observational signature, the total duration of inflationary epoch should not be much longer 
than the observable e-foldings A'^obs. • This tuning represents a highly non-trivial constraint 
on inflationary physics. 

It would be interesting to extend our results to generic multiple field models where the 
long-wavelength effects could be even more pronounced, as we have briefly discussed above. 
Another interesting line of future work would be to consider non-Gaussianities generated by 
subhorizon physics which generate non-trivial deviations from the locai form. Developing 
a solid understanding of the long-wavelength effects would indeed be of key importance as 
it appears to play an essential role in translating the new observational data into a deeper 
understanding of fundamental high-energy physics. 



Note added Whilst this paper was being completed, a related work appeared on the 
arxiv by Nelson and Sbanderà |26]. This work is conceptually similar to ours, since it also 
studies the relation between our observable correlators to those in a much larger patch 
in single-source scenarios and focusing on the locai model. However their focus is on the 
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structure and hierarchy of the non-Gaussian moments, while we focused on more concrete 
calculations, such as explicitly finding how much /nl varies under the influence of 5nl- 
This paper and |26| are therefore complementary. 
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